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Time-periodic perturbation can be used to modify the transport properties of the surface states of 
topological insulators, specifically their chiral tunneling property. Using the scattering matrix method, we 
study the tunneling transmission of the surface states of a topological insulator under the influence of a 
time-dependent potential and finite gate bias voltage. It is found that perfect transmission is obtained for 
electrons which are injected normally into the time-periodic potential region in the absence of any bias 
voltage. However, this signature of Klein tunneling is destroyed when a bias voltage is applied, with the 
transmission probability of normally incident electrons decreasing with increasing gate bias voltage. 
Likewise, the overall conductance of the system decreases significantly when a gate bias voltage is applied. 
The characteristic left-handed helicity of the transmitted spin polarization is also broken by the finite gate 
bias voltage. In addition, the time-dependent potential modifies the large-angle transmission profile, which 
exhibits an oscillatory or resonance-like behavior. Finally, time-dependent transport modes (with 
oscillating potential in the THz frequency) can result in enhanced overall conductance, irrespective of the 
presence or absence of the gate bias voltage. 



Topological insulators (TI) form a new class of quantum matter possessing insulating bulk and metallic edge 
or surface states''^. The surface states of TIs have an odd number of massless Dirac cones', which are 
protected against time-reversal-invariant perturbations such as non-magnetic impurities, defects, and 
reconstruction. This promising characteristic motivates significant interest in the fundamental physics of trans- 
port in TIs"* Due to the robustness of the surface states against time-reversal symmetric perturbations, TIs are 
expected to be good candidates for the design of spintronics devices'". TI is also a good candidate for fault- tolerant 
quantum computing" because of the formation of Majorana bound state induced by the superconducting 
proximity effect'^. 

To realize the practical Tl-based devices, one major task is to effectively control and adjust the transport 
properties of the surface states. Correspondingly, some interesting effects, such as finite-size effects " ferromag- 
netic coupling'^'"", and quantum confinement" '" have been utilized to achieve this goal. In semiconductors one 
can confine electrons quantum mechanically using electrostatic barriers, e.g. a quantum well. However, such 
electrostatic barriers cannot provide electron confinement in graphene and TI systems, due to chiral tunneling of 
Dirac fermions, a phenomenon which has been extensively studied recently""^^. The chirality of electron's 
wavefunction in graphene gives rise to perfect quantum tunneling for electrons incident in the normal direction 
of a potential barrier^"". The chiral nature of Dirac fermions also suppresses the backscattering of quasiparticles 
leading to high charge carrier mobOity, which can be utilized in some novel nanodevices^'''^'. Similarly, the surface 
states of TIs can also mimic Dirac fermions in quantum electrodynamics (QED). Naturally, it is meaningful to 
study the chiral tunneling of quasiparticles when a potential barrier is applied to the topological surface states. In 
recent years, much attention is paid to optical response of the topological surface states^"". Applying a time- 
dependent perturbation can provide an alternative and efficient way to control spectrum and transport prop- 
erties, e.g., inducing a topological state in a semiconductor quantum well" and a conventional insulator'^, 
resulting in Floquet-Bloch bands due to the hybridization of bulk states and topological surface states'"". 
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Figure 1 | Schematic of the conducting channel with time-periodic 
potential comprising of topological surface states. 

However, relatively few studies have been conducted on the chiral 
tunneling of surface states under the influence of a time-dependent 
potential, or under nonequUibrium condition, i.e., in the presence of 
a gate bias voltage. 

In this work, we consider a conducting channel consisting of 
topological surface states influenced by a time-periodic potential, 
shown schematically in Fig. 1. It is found that in general, a time- 
dependent potential modifies the large-angle transmission of the 
device, which exhibits oscillatory or resonance-like behaviors. On 
the other hand, the presence of a gate bias voltage can significantly 
affect the transmission probability of normally-incident electrons, 
and destroys the perfect transmission for normally- incident elec- 
trons, which is a signature of the Klein tunneling effect^''. The char- 
acteristic left-handed helicity of the transmitted spin polarization is 
also destroyed by the finite gate bias voltage. However, the zero-angle 
spin polarization is robust against the time-dependent potential and 
the gate bias voltage. Finally, by applying the scattering matrix 
method, we numerically study the amplitude and frequency depend- 
ence on the transmission probability, and describe the time-depend- 
ent response of the systems conductance in the presence or absence 
of a gate bias voltage. 

This paper is organized as follows. In section Model and formula- 
tion, we consider the Dirac equation governing the topological 
surface states of the system in the presence of the time-periodic 
potential in the central region and a gate bias voltage across the left 
and right portions of the conducting channel. By solving the Dirac 
equation, we obtain the analytic solution of the time-dependent 
wavefunction. We employ the scattering matrix method to calculate 
the transmission probability and the device conductance. In section 
Numerical results and discussion, we analyze the dependence of the 
transmission probability on the electron injection-angle and the fre- 
quency dependence of both the transmission probability and the 
overall conductance. We summarize the main results in section 
Summary. 

Model and formulation 

In this section, we investigate the transport property of a conducting 
channel proposed of topological surface states, as illustrated in Fig. 1, 
modulated by a time-dependent potential, which arises from a small 
ac signal V{t) applied to the gate voltage in the central region. 

At low energies, the surface states can be described by the effective 
HamUtonian' 



H=Vp[pya^-p^(Ty), 



(1) 



where Vp is the velocity of the surface states, vp ^ 5 x lO'm/s for the 
typical topological insulator 612863. The wavefunction of the time- 
dependent system satisfies the following Dirac equation 



- ihvp 



0 



il/^-li{x)il/f = ih8tilJi, 



where i = 1, 2, 3 are the notations of three different regions, 
xj/. = is the two-component wavefunction with f, J, being 

spin indexes. The potential function /((x) is written as 



Hi^, x<-L/2, 
li{x) = l —UQCOs{ujt), —L/2<x<L/2, 
l^ij^, x>L/2. 



(3) 



This potential profile of the conducting channel can be envisaged by 
applying top-gate and back-gate voltages^'""*", where Ha^L,R describe 
the relative shift of Dirac points of region I and III with respect to 
region //, V = (/i^. ^ /ijj)/e is defined as the gate bias voltage, and Uq 
and oj are the amplitude and frequency of the time-dependent 
potential. 

Firstly, we try to find the solution of Eq. (2) for the central region. 
The wavefunction can be constructed as*' 



GO 



(4) 



where i/'2o('') is also a two-component wavefunction, whose explicit 
form is given below, and /;„ refers to the mth-order Bessel function. 
Substituting Eq. (4) into Eq. (2), we can obtain the equation satisfied 
by t/'2o(r) 



- ihvp 



dy-id^ 



0 



1/^20 (r)=£i/'2oW 



(5) 



Thus the solution of the two-component wavefunction can be 
expressed as 



^20,+ ^ 



v/2 



ye 



1 



(6) 



where k=yjk^ + k^ is the amplitude of the wave vector, 

0 = !LTcsm{hvpky/\e\), kx(e,) — \£\cosO/{hvp), and )'(e) = sgn(£) = 
±1 corresponding to electronlike and holelike cases, respectively. 
Correspondingly, the wavefunctions i/'i3(r, i), in the regions j^j > 
L/2, obey similar formulas to Eq. (5), and the solution of which can be 
written as 



1 

V2 



where q=yqx~^'iy i^ ^1^° '^e amplitude of the wave vector, 
Oo.^. = arcsm{hvpqy/\e + ^iJ), qx,a{ii) = \e + ^iJcosOo^a/{J^vp), and 

yo.a(«) = sgn(e + ^la) = ±1. 

Thus for a charge particle with an energy £ which is injected from 
region I into the center region, the wavefunctions in the three regions 
can be written as 



yo.il 



1 

-|-mfiffl)e'("i'"-°i" 
1 



+ > . '■m(£) 

m— — CO 



(8) 



x< 



-L/2, 
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+ B(m) I |e 



-L/2<x<L/2, 



5Z 



>'o_j;(e + mfifo) £'(2 -"(»".«).«) 



x>L/2. 



(10) 



Imposing the boundary conditions for the wavefunctions i/',(r, t), 
which must be continuous at the interfaces of x = ±L/2, we can 
obtain the equations relating A{m') and B(m') utilizing the ortho- 
gonality of {e^""'"'}. After some algebra, we obtain the coefficients 
A{m') and B{m'), and substitute them into the following equation to 
calculate the transmission coefficients 



tm e =e 



- iqxR{i^ + f^^('j)L/2 



J2 [Aimy'' 



{E + m'tifo)L/2 



(11) 



After these quantities are obtained, we can study the quantum trans- 
port property of the considered system using the scattering matrix 
method*^ Correspondingly, we can obtain the total transmission 
probability 



,„ HxL\<-) 



(12) 



Consider an electron incident from region / with a fixed Fermi energy 
Ep. The electron can be injected into the central region at various 
incidence angle 0, associated with different y-component wavevector 
q_yo- We can calculate the transmission probability for a fixed q^o or ^' 
i.e., T{0, Ep). The overall zero-temperature conductance G can be 
calculated by integration of Eq. (12) over they-component wavevec- 
tor"'''^, which can be written as 



2e^ W 
h 2nh 



dpyT{py)=Go dOcosOT{0,Ep), (13) 

J J -k/2 



where the factor 2 accounts for the twofold spin degeneracy, W is the 
transverse dimension of the interface and Go = e'kpW/{nh) is the 
unit of conductance. 

It is well known that the electron spin is locked in the surface plane 
for a purely 2D system*"-'^"". Naturally, it is interesting to study 
whether this spin-locking characteristic can be maintained in the 
presence of the time-periodic potential and the gate bias voltage. 
Based on the wavefunctions in Eqs. (8) to (10), we can calculate 
the spin polarization averaged over spin-space, i.e., Pf(f) = 
H/2(i/'|(T,,|i/'), with € = x,y, z. The time average of the spin polariza- 

tion is then defined as P( = oj I dtV({t) j2n. 



Numerical results and discussion 

In this section, we numerically study the transport property of the 
topological surface states modulated by a time-dependent potential 



for varying gate bias potential. We first analyze the total transmission 
T{0, Ep) as a function of the electron injection-angle 0 for different 
amplitudes and frequencies of the time-dependent potential, as 
shown in Fig. 2. In our numerical calculations, the following para- 
meter values are assumed: Fermi energy of Ep = 100 meV, Fermi 
velocity of = 5 X 10^ m/s, and length of central region (where the 
time-dependent potential is applied) of L = 60 nm, which provides a 
sufficiently thick barrier for electrons with a Fermi wave length Ap ~ 
20 nm. Under zero shift of the Dirac points in region I and III, i.e., 
A'l = /^j; = 0, virtually perfect transmission {T ~ 1) of electrons 
through the time-dependent central region occurs for incidence 
angles \ 0\ < 70°, (see Fig. 2 (a)). Such perfect transmission of normally 
incident electrons is a well-known feature of Klein tunneling of Dirac 
fermions^' One can understand the perfect transmission as arising 
from the perfect matching between the electron wavefunctions of the 
three different regions (see Eq. (6), (7)). In the presence of a time- 
dependent potential, the large-angle transmission of electron is 
modulated, and the modulation becomes more obvious at larger 
amplitude, e.g., Ug = 3hoj represented by the black curve. This phe- 
nomenon can be ascribed to the time-dependent potential, as there 
exists no other potential barrier over the entire surface of the device. 
In Fig. 2 (b), we applied gate voltage of the same sign to the left and 
right parts of the conducting channel, i.e., f ip = /.i^ = 50 meV. In this 
case, the large-angle transmission is suppressed, and electron trans- 
mission is mainly confined to a smaller range of incident angles of 
|0| <40°. However, perfect transmission is still observed at normal 
incidence, i.e., at f = 0. The effect of the time-dependent potential on 
the transmission probability is relatively minor and confined to some 
oscillations at 0 close to 40". When the applied gate voltage is 
increased, i.e., to /.ip = fiR = 200 meV, the transmission is further 
confined to a smaller range of 0, as shown in Fig. 2(c). Perfect trans- 
mission, however, still persists at normal incidence. 

It is interesting to note perfect transmission of normally incident 
electrons no longer occur when a finite gate bias voltage is applied. 
Even a small bias voltage, e.g., y = 10 meV, can result in the dis- 
appearance of the perfect normal tunneling (see Fig. 2(d).) When an 
electron is injected normally from region I, and transmits through 
the central region into the region III, the asymmetry in the gate 
potentials of region / and III destroys the perfect matching of the 
electron wavefunctions in the three regions, resulting in a finite 
reflection and hence imperfect transmission process. This phenom- 
enon is a striking result in this paper. In addition, the transmission 
profile is slightly modulated by the time-dependent potential, espe- 
cially near the threshold values of 0 where there is a sharp drop in the 
transmission probability. The effect of the time-dependent potential 
becomes weaker with increasing gate bias voltage. As shown in 
Fig. 2(e), when V = 100 meV, the modulation due to the time- 
dependent potential is hardly visible. Next, we study the effect of 
the frequency of the time-dependent potential on the transmission 
profile. The transmission T{0, Ep) is plotted for different frequencies 
with a fixed value of UQ/hco= 1.5 (see Fig. 2(f)). It is observed that 
high-frequency potential induces a stronger modulation of the large- 
angle transmission, as reflected by the larger oscillations in the trans- 
mission profile. 

To further analyze the effect of the time-dependent potential, we 
investigate the transmission probability as a function of the fre- 
quency for different amplitudes with a fixed incident angle, as shown 
in Fig. 3(a). Under equal gate voltage and no bias condition of fip = 
/.Ir = 50 meV, the transmission shows oscillations with frequency, 
with the oscillatory behavior becoming more pronounced with 
increasing amplitude of the time-dependent potential. Next, we con- 
sider the asymmetrical case, f-ip = ^^Ir = 50 meV, as shown in 
Fig. 3(b). As in the previous case, the time-dependent potential 
induces an oscillatory behavior in the transmission with respect to 
frequency. However, the transmission probability is significantly 
smaller than that in the equal voltage case, and drops to almost zero 
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Figure 2 | The total transmission T(0, Ep) plotted as a function of the electron injection-angle 6, for (a) = p.K = 0, (b) 50 meV, and (c) 200 meV, 
for (d) /(i = = 5 meV, (e) = — /ij; = 50 meV. Green, blue, red and black curves correspond to Uq = 0, 0.5, 1.5, 3 fico with to = 3.05 THz. 
Especially, in the subgraph (f), = 5 meV, and Green, blue, red and black curves correspond to co = 3.05, 6.1, 9.15, 12.2 THz with i<o = 1.5 fioj. 

The Fermi energy is chosen to be Ef = 100 meV and the Fermi velocity is Vp = 5 X 10^ m/s. 
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at certain range of frequencies, e.g., 21a)o <oJ<25a)o. We also find 
that the transmission probability increases with increasing amplitude 
of the time-dependent potential. By contrast, in the equal-voltage 
case, the reverse trend occurs, i.e., the transmission probability is 
lowered as the potential amplitude Uq is increased for the incident 
angle 0 = 30°. For smaller incident angles oi\0\ < n/lO, the trend of 
increasing T with increasing Mq remains. Since the bulk of the trans- 
mission occurs with electrons having small incident angles of |0| < 
71/ 10, the overall conductance increases with increasing Mq in the 
equal- voltage case. 

In Fig. 3, we analyze the transport behaviors only for a specific 
incident angle. Next, we consider the effect of the time-dependent 
potential on the overall zero-temperature conductance G, which 
comprises of contribution of all electrons incident within the angular 
range of —K/2<0<n/2. In Fig. 4, G is plotted as a function of 
frequency co for different potential amplitudes Uo. One can see that 
the conductance G is significantly higher in the equal-voltage case 
(Fig. 4(a)), compared to the asymmetrical case (Fig. 4(b)). However, 
for both cases, the conductance G increases with the potential ampli- 
tude Mq. Thus, one can surmise that time-dependent transport modes 
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Figure 5 | Tlie time-averaged spin polarization P plotted as a function of 
tlie electron injection-angle 0 for m = 3.05 THz and uq = Ihco. Other 
parameters are the same as those in Fig. 2. 



(with oscillating potential in the THz frequency) can result in 
enhanced overall conductance, irrespective of the presence or 
absence of gate bias voltage. 

To analyze the spin-locking characteristic of the spin polarization, 
we investigate the time-averaged spin polarization P = P/ex + Pyey + 
P^Cz at X = L/2, as a function of the electron injection-angle 0 for 
different gate bias voltages, as shown in Fig. 5. The length of the 
arrows denotes the amplitude of the spin polarization, while a dot 
represents zero spin polarization. In the numerical calculation, the 
Fermi energy is set to Ep = 100 meV, the frequency and the ampli- 
tude are frxed to be co = 3.05 THz and Mq = 2fioj. It is found that the 
z-component spin polarization is zero, P^ = 0, i.e., the spin polariza- 
tion is still locked to the x-y plane. When there is no shift of Dirac 
points in region I and J//, i.e., ^ii = jij^ = 0, P has a left-handed 
helicity and is always perpendicular to the wavevector qo of the 
incident electrons. For this case, the spin polarization exists over a 
large angular range, 0 < 70", which is consistent with the large-angle 
perfect transmission in Fig. 2. Under the equal-voltage case, e.g., = 
= 50, 200 meV, P is still perpendicular to the wavevector qo, but 
the angular range of the spin polarization decreases with increasing 
gate potential. Thus, one can adjust the spin polarization and confine 
it to the small-angle region by changing the gate voltages of the left 
and right portions of the conducting channel. Interestingly, we also 
observe that the characteristic left-handed helicity of the transmitted 
electrons is also broken in the presence of a finite gate bias. The 
deviation from the helical spin orientation becomes more pro- 
nounced as the bias voltage is increased (as can be observed by the 
green arrows corresponding to the largest applied bias). In addition, 
we found that the deviation from the left-handed helicity is greater 
for large-angle transmission. For the zero-angle transmission, the 
spin polarization is robust with respect to the applied bias voltage 
and time-dependent potential. In terms of the definition of the spin 
polarization, when 0 = 0°, we find that Px/y{t) «: J2m,m' C(^)*m'(£') 
fm[}'o,R{'^ + fn'fK^)+}'Q,R{£ + fnhco)], where is an exponential 
function including m and m' . Basically, the condition 

n{£ + m'hco) —y'Q h{e + mhoj) =0 can be satisfied for the small 
values of m and m' in the above five cases, while the transmission 
coefficients will become zero for larger values of m. Thus for the 
special case oi 0 = 0°, Pjdt) becomes zero and this results in the 
robustness of the spin polarization. 

Summary 

In the above, we present a theoretical study of the influence of time- 
dependent potential and gate bias voltage on the chiral tunneling for 
topological surface states. Utilizing the scattering matrix method, we 
numerical analyze the effect of the time-dependent potential on the 
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transmission probability, the conductance and the spin polarization. 
We find that, under the equal- voltage case (where the gate potentials 
of region I and III are identical), perfect transmission occurs for 
electrons which are injected normally to the interfaces. The spin 
polarization retains its left-handed helicity, but the angle range of 
the spin polarization decrease with increasing gate voltage. However, 
in the presence of an asymmetrical gate voltage two signature char- 
acteristics of Klein tunneling no longer hold. Firstly, normally incid- 
ent electrons no longer experience perfect transmission-the 
transmission probability of such electrons diminishes with increas- 
ing gate bias voltage. Secondly, the left-handed helicity of the spin 
polarization of transmitted electrons is also broken by the finite gate 
bias. As for the time-dependent potential, its influence is prominent 
only for large-angle transmission, where it induces oscillatory or 
resonance-like behavior in the transmission. The overall conduc- 
tance of the device is significantly reduced upon application of a gate 
bias voltage. Finally, the presence of time- dependent transport 
modes which give rise to time- dependent potential, would enhance 
the device conductance, irrespective of the presence or absence of 
gate bias voltage. 
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